NON-PERTURBATIVE MASS AND CHARGE 
RENORMALIZATION IN RELATIVISTIC NO-PHOTON 
QUANTUM ELECTRODYNAMICS 



CHRISTIAN HAINZL AND HEINZ SIEDENTOP 



Abstract. Starting from a formal Hamiltonian as found in the physics lit- 
erature - omitting photons - we define a renormalized Hamiltonian through 
charge and mass renormalization. We show that the restriction to the one- 
electron subspace is well-defined. Our construction is non-perturbative and 
does not use a cut-off. 

The Hamiltonian is relevant for the description of the Lamb shift in muonic 
atoms. 



1. Introduction 

According to Dirac's hole theory the vacuum consists of electrons which occupy 
the negative energy states of the free Dirac operator (Dirac sea). Dirac postulated 
that their charge is not measurable. However, if one introduces an external elec- 
tric field, e.g., the field of a nucleus, these electrons should rearrange, occupying 
the negative energy states of the Dirac operator with the external electric field. 
Physically speaking, the nucleus polarizes the vaccuum. (This rearrangement may 
be interpreted as the creation of virtual electron-positron pairs when expressed in 
terms of the free Dirac operator.) In other words, the vacuum is polarized. Dirac 
indicates that these polarization effects result in a logarithmically divergent charge 
density, which cannot be neglected. As a solution, he suggested that a momentum 
cut-off must be introduced, since he expected that the Dirac equation would fail for 
energies higher than 137mc 2 . In he changed this train of thought and suggested 
that the infinities occurring should be absorbed by a procedure which is now called 
charge renormalization. A similar step was independently undertaken by Furry and 
Oppenheimer in jH| who circumvented the hole theory by introducing annihilation 
and creation operators. 

Heisenberg clarified Dirac's picture and generalized his approach extracting 
the physically relevant terms by subtraction of an unambiguous infinite constant, 
at least to first order in a. Serber [23 and Uehling gave detailed calculations 
(in first order of a). Uehling demonstrated that the vacuum polarization alters the 
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Coulomb potential of a charged particle resulting in the electron being slightly more 
bound in the s-states (angular momentum 0) of hydrogenic atoms. Later Weisskopf 
gave a thorough discussion of the physics involved in charge renormalization. 

After the experiments of Lamb and Retherford j^j in 1947, which gave a much 
higher discrepancy concerning the hyperfine structure of hydrogen, in addition to 
a different sign, than Uehling's calculation showed, and the first explanation by 
Bethe pQ, the insight into quantum electrodynamics (QED) changed and the inter- 
action with the radiation field turned out to be the dominating part in describing 
the splitting of the energy levels of hydrogenic atoms beyond the Dirac equation. 
Similar to vacuum polarization, which was now treated together with the radiative 
corrections, the photon interaction caused fundamental problems such as infinities, 
which were "removed" - at least in first order of a - by mass renormalization 
by Tomonaga, Schwinger, and Feynman. Eventually, Dyson "succeeded" with the 
renormalization program to every order in a. Since then, QED has proven to be of 
extraordinary predictive power. (We refer the reader interested in more historical 
details to Schweber \'24\.) 

But despite the predictive power of quantum electrodynamics, the description 
in terms of perturbation theory causes great uneasiness among mathematicians; a 
mathematically consistent formulation of QED is still unknown; in fact Dyson |S] 
indicated that the perturbation theory is divergent. A self-adjoint Hamiltonian for 
QED is not known. 

In the present paper we address a particular kind of singularities arising in QED, 
namely those stemming from the vacuum polarization. As opposed to the prevalent 
physics literature we will not use any Feynman diagram but will rather construct 
a Hamiltonian (in Coulomb gauge). This we have in common with the above cited 
early works in the field. However, the fact that we start from a formal Hamiltonian 
and renormalize it non-perturbatively distinguishes us from those authors. 

Altough our approach is rigorous, the resulting renormalization is far from being 
of academic interest only. In fact, the restriction D rcn of the fully renormalized 
Hamiltonian H to the one-particle electron sector accounts already for a precise 
description of the low energy levels of /i-mesonic atoms where the vacuum polariza- 
tion effect dominates the radiative corrections by far, since the Bohr orbits traverse 
the support of the polarization potential in this case. (See, e.g., Peterman and 
Yamaguchi Glauber et al. Milonni ^J, Weinberg POJ, and Greiner et al 

HBO 



2. Model 

In relativistic QED the quantized electron-positron field \P(x), which is an oper- 
ator valued spinor, is written formally as 

(1) = a(x) +b*(x) 

where a(x) annihilates an electron at x and b* (x) creates a positron at x. (We use 
the notation that x = (x, <r) 6 T = M 3 x {1,2,3,4}, where a is the spin index and 
/ dx denotes integration over M 3 and a summation over a.) The underlying Hilbcrt 
space is given by S) = L 2 (T). 
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The definition of a one-electron, respectively one-positron, state will correspond 
to the positive, respectively negative, energy solutions of the Dirac operator 

(2) D v = a ■ iv + (3-a<p 

i 

in which a,/3 denote the 4x4 Dirac matrices. The constant a is a positive real 
number, the Sommerfeld fine structure constant which is approximately f/f37. (We 
have picked units in which the electron mass is equal to one.) 

We will not assume that the nucleus is a point particle; we rather associate with 
it a density n 6 L 1 (M) whose integral gives the atomic Z number of the atom under 
consideration. For technical convenience we assume n to be a spherically symmetric 
Schwartz function whose Fourier transform has compact support. 

We remark: it is an experimental fact that the nucleus is not a point particle but 
an extended object. In fact the numerical calculations of the Lamb shift depend on 
the size of the support, which actually limits the accuracy of numerical value of the 
calculation of the Lamb shift because of the experimental uncertainty of the radius 
of the nucleus (Weinberg 30 , p. 593). (A point nucleus leads also to mathematical 
difficulties, since the renormalized potential is more singular than the Coulomb 
potential, i.e., it could not be controlled by the kinetic energy (see Uehling and 
Subsection 13 . 51 of this paper). 

The electric potential of the nucleus is given as 

(3) v = \-\~ 1 *n. 

An application of the Young inequality shows that the nuclear potential 

(4) ip Gi 3+e (M 3 ) ni°°(m 3 ) 

under our assumption on the nuclear density n for any positive e. Moreover, by 
Newton's theorem 

(5) < p(x) < Z/|x|. 

For completeness we note the following fact whose proof is obvious: 

Lemma 1. Fix aZ € [0,1) and assume n to be a non-negative Schwartz function 
with J n = Z and tp = | • | _1 * n. Then, D v is selfadjoint with domain H 1 ^), 
i.e., the same domain as the free Dirac operator Dq. Moreover, if n is spherically 
symmetric, D v has no negative discrete spectrum and its k-th eigenvalues is bounded 
from below by the k-th eigenvalue of the Coulomb Dirac operator D z l\'\. 

We remark that Lemma ^ implies that the lowest eigenvalue is positive. 

Now, we can specify the electron and positron state spaces $j+ and CS}f_ respec- 
tively: the orthogonal projection on Sj+ and CSjf. are defined as Pf :— X[o,oo) (f„. m ) 
and Pf_ := 1 — P+. The (anti-unitary) charge conjugation operator is given on F) by 

(Cip)(x) = i{3a2ip(x). In momentum space it acts as (Cip)(p) — i(3c<2ip(—p)- (Here 
and in the following we follow the notation of Thaller [2HI; see also |12|.) 
Note that the Hilbert space can be written as the orthogonal sum 

(6) f)=S%.®Sst. 

Correspondingly a*(f) creates an electron in the state P+f, whereas b*(g) creates 
a positron with wave function CPf_g. Note that the definition of the operator a and 
b depends explicitly on the choice of the potential tp. 
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The Hamiltonian for the non-interacting electron-positron field is given by 

(7) »-./* ..www 

where : : denotes normal ordering, i.e., anti- commuting of all creation operators to 
the left of all annihilations operators ignoring the anti-commutators. 

Note, that for our renormalization procedure the choice of the electron-positron 
subspaces as the positive and negative eigenspaces of D v is crucial, in fact it is a 
choice already proposed by Dirac 

The creation and annihilation operators fulfill the canonical anti-commutation 
relations 

{«(/), a(ff)} - K(/), «*(<?)} = {a(f),b(g)} = {«*(/), b*(g)} 
{S) ={a*(f),b(9)} = Mf),b*(9)} = 0, 

and 

(9) {«(/), a* (g)} = (f,P£g), {&*(/),%)} = {f,P?g). 

Formally this is equivalent to 

{a(x),a(y)} = {a*(x), a*(y)} = {a(x),b(v)} = W (x) , b* (y)} 
= {a*(x),b(y)} = {a(x),b*(y)} = 0, 

and 

(11) {a(x), a*(y)} = P*(x, y), {b* (x), b(y)} = P?(x, y) 

where Pf (x, y), Pf.(x, y) are the integral kernels of the projectors Pf, Pf_. 

If there is no external potential there should be no polarization effects present. 
It is therefore expected that the difference Q v of one-particle density matrices of 
the perturbed and unperturbed vacua 

(12) Q v := Pi - P { 1 = -P* + P°_ 

plays a central role in defining the renormalizcd Hamiltonian. Using Cauchy's 
formula we can express the Q v in terms of the respective resolvents (Kato |14| . 
Section VI, 5, Lemma 5.6) 

1 f°° ( 1 1 

(13) Q v = — dq' 



(10) 



2-kJ_ 00 yD^ + iri D a + iri / 
The difference of Q v and of the first order resolvent expansion, i.e., 

/°° 1 1 

will contribute to the renormalized operator as follows: one interprets its spin 
summed diagonal as density. The corresponding electric potential should be added 
to the one particle operator. To avoid any unnecessary difficulties defining the 
operator we split (|14(l again in three summands motivated by iterating the resolvent 
equation: 

(15) a 2 Q 2 + a 3 Q 3 + a 4 Q 4 , 
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where the indices 2,3,4 indicate the number of </?'s in the expression, i.e., 



Q 2 := 



1 



(16) 



1 



dr\ 
dr\ 



1 



1 



1 



D° + irj D" + if) D° + irj' 
1 1 1 



-ip- 
-ip- 
-ip- 



-ip- 
-p- 
-p- 



i 



D° + hf D" + ir\ r D° + irj^D + irj' 

1111 

-<p- 



1 



-tp- 



D° + irj D° + irj D<p + irj D° + irj D Q + irj 

We can immediately remark that the density corresponding to Q 2 vanishes: the 
terms linear in the Dirac matrices vanish after summation over a, since the Dirac 
matrices are traceless; the remaining terms are odd in 77 and vanish after integration 
over 77. Consequently, we can disregard this term in defining the operator. 

We now define the density 

4 



(17) 

where 
(18) 



Ps(x) := (2 



ny 3 f dp f d q y2e^-^Q 3 (p,a;q,a) 

Jm. 3 Js. 3 rT=1 



Mp,q) = 7T 



dr] / dpi / dp 2 (D 



irj) 1 o 



^(P - Pi) ° (-Dpi 



o (p(p 1 - p 2 ) o (D P2 + irj) 1 o p(p 2 - q) o (_D q 



ir,)- 1 

-in)- 1 



a • r + (3. The corresponding electric potential is 

-1 



with D, 

(19) P3 := P 3 * 
The quadratic form defining P4 is given by 

(20) W,Pi4>) :=tr( X Q 4 ) 

where x(x) := / dy\\j){y)\ 2 /\y - x|. 

It will be useful to introduce the function C 



(21) 



1 f 1 

C(k) = -k 2 / dx(l - x 2 ) log[l + k 2 (l 

2 Jo 



! )/4] 



h 2 i L 4", Vl + 4/k 2 + l 



4 

k 2 " 



k 2 ~~° ^/l + 4/k 2 - 1 

as already done by Serber [SSj and Uehling and later by Pauli and Rose 20 , 
Jauch and Rohrlich J3|> Schwinger [25], and Klaus and Scharf 16 . The vacuum 
polarization potential U, also known as Uehling potential, is defined via its Fourier 
transform 

,C(k) 



(22) 
T 

(23) 
whe 

(24) 



U(k) = 0(k) 

The renormalized one-particle operator is 
1 



7r|k| ; 



D rcn := a. ■ tV + - onp - a 2 U + a 2 X - a 4 P 3 - a 5 P 4 



where X is the renormalized operator with integral kernel 

<yp{x, v ) 



aX(x,y) := 
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To introduce the operator D ren might appear unmotivated at this point. How- 
ever, it has a solid physical motivation: it emerges through mass and charge renor- 
malization from the canonical formal text book Hamiltonian (see, e.g., Milonni 19 , 
p. 385, Formula (11.25)) for the interaction of electrons when there are no photons 
present. In turn the charge and mass renormalization originates in three physical 
principles Wl, W2, and W3 as we will explain in Sectional 

Moreover, and this is our main mathematical result, the operator P^D rcn on f)+ 
turns out to be well defined and self-adjoint on the same form domain as the free 
Dirac operator: 

Theorem 1. Assume aZ G [0,1) and ip — n * \ ■ | _1 with spherically symmetric 
Schwartz function with compact support in Fourier space. Then the quadratic forms 
(ip,Uip), (ip,Xi/j), (ip,P 3 ip), (%p,P^ip) on P^(_ff 1 (r)) are relatively form bounded 
with respect to (ip, l-D ^) with form bound zero. 

This has the following consequence: 

Corollary 1. There exists a unique self-adjoint operator D + fulfilling (ip, D + ip) = 
(ip, (D v + a 2 X - a 2 U - a i P 3 - a 5 Pi)ip) for all ip G Pf (if 1 ^)) with form domain 
P^_(H 1 ^ 2 {T)). Furthermore it is bounded from below. 

Proof. Define c := \/2 and M := V^IMU- Then, obviously < (c 2 - 2)(L> ) 2 - 
2tp 2 + M 2 implying (D ) 2 < c 2 (D 2 - 2$((pD ) + (p 2 ) + M 2 + 2Mc\D% Rewriting 
this and taking the square root - which is an operator montone function - yields 

(25) \D°\ < c\D v \ + M. 

This means that infinitesimal form boundedness with respect to \D°\ implies also in- 
finitesimal form boundedness with respect to \D V \. Moreover, we have (ip, \D v \ip) = 
(ip,D^ip) for ip G P^H 1 ^)). 

Thus, according to the KLMN theorem (Reed and Simon [231 , Theorem X.17) 
there exists a unique self-adjoint operator D + on io^ whose form domain is the 
form domain of D v (which - since <p is bounded - is the form domain of P+D°, 
i.e., P+if 1/,2 (r)). Moreover, the quadratic form of D + fulfills 

(26) (ip, D + ip) = (ip, (D* + a 2 X - a 2 U - a 4 P 3 - a 5 P 4 )^). 

□ 

As already mentioned in the introduction, the Hamiltonian (|23[) can be used to 
describe /j,-mesonic atoms where the interaction with the photon field is negligible 
as indicated experimentally by Peterman and Yamaguchi |21) and theoretically by 
Glauber et al [5]. 

3. Physical "Derivation" of the Renormalized Hamiltonian 

We start with the formal expression for the interaction of electrons when no 
photons are present (Kroll and Lamb ^I|> French and Weisskopf 0) found also in 
text books (see, e.g., Milonni |19j): 

(27) y fw = iUf dy ^E}^m_ 

2 7 J pc-y| 

The Hamiltonian describing our system is formally given by 

(28) H ur = B£ r + aW UI 
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where 

(29) nt^/d.vw^.w 

It is well known that this expression contains several singular terms. In particular 
it does not even contain any normal ordering. The remaining part of this section 
can be viewed as manipulating on it a physical allowed way and transforming it to 
a physically equivalent expression that is mathematically meaningful, namely the 
renormalized Hamiltonian. We emphasize that none of the steps taken is mathe- 
matically justified, i.e., the eventual justification of the renormalized Hamiltonian 
is its successful predictive power. 

We use three guiding principles to transform expressions for the energy into other 
physically equivalent ones as formulated and justified by Weisskopf P- 6: "The 
following three properties of the vacuum electrons are assumed to be irrelevant: 

Wl: The energy of the vacuum electrons in field free space. 

W2: The charge and current density of the vacuum electrons in field free 
space. 

W3: A field independent electric and magnetic polarizability that is constant 
in space and time. " 

Similar procedures have been suggested by Heisenberg JT] , French and Weisskopf 
0, Kroll and Lamb [T7j, and Dyson g]. 

Exploiting the canonical anti-commutation relations (jlOJI we can rewrite (12811 . 
For the one-particle part we have 



(30) = W + J dx(D v P^)(x, x). 

The last summand is a - although infinite - constant which we drop, since it does 
not influence energy differences. For the two-particle part we get 



(31) 



f J , T ,„ > Tr , , P?(x,y) - P*(x,y) 



^Jdxjdy: $*(x)V(y) 



P v (v v) 

dx I dy: • ~ Ky,y> 



|x-y| 

If f J P%(x,y)P*(x,y) If f P*{x,x)P*{y,y) 

+ 2j dX J dV ^^\ + 2 J dX J dy ^^V\ ■ 

The last two terms are again constants which we drop. The first term is the 
normal ordered two-particle interaction which has finite expectation in states of 
finite kinetic energy. We will denote it by 

(32) W=l[ d x[ d y^ X ^*^-- 



ix-yi 

The remaining two other terms are one-particle operators of particular interest to 
us. 

Both terms, the classical electrostatic interaction energy of the electron with the 
polarized Dirac sea called the "non-exchange energy" 

(33) Pur :=-JdxJdy. **(x)*(x) : ^z^f , 
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and the exchange energy 

(34) X m . :=Ud X [dy: ' " 



2 J J w ya ' |x-y| 

are not well defined. (For curiosity we remark that the latter is logarithmically 
divergent in A, if one introduces a cut-off by ^>{x) = J| p i< A ^(p, er)e _ipx dp). 

To renormalize the exchange energy we introduce the operators P^_,P®_ which 
are the projectors on the positive and negative subspace of the free Dirac operator 
D°. (Note that we can interpret P° as the one-particle density matrix of the free 
Dirac sea.) 

3.1. The Renormalization of the Exchange Energy. To renormalize X ur we 
subtract the exchange interaction energy of the electron with the free Dirac sea 
using Principle Wl, i.e., 

1 [j. f ,...«.*,^„.,. [(PZ-PZ)-(P*-P°)](x,y) 



aX :=~ f dx f dy : ^*(x)^(y) 

(35) 



dx / dy : : X(x,y) 



with X as defined in (|2"H) . 

(In physical language this subtraction of an undefined operator - known as 
"counter term" - is called "mass renormalization" . We refer to French and Weis- 
skopf 7 , Equation (30), for the motivation of this terminology.) 

From now on, we will assume that the external potential if is so weak that all 
there are only positive eigenvalues in the gap of D v . Then, 

1 1 r°° 1 

(36) Pi = - + — / dn 

as well as 

1 1 

2 + 2^ y_ OQ ' D° + irj 
(Kato Section VI, 5, Lemma 5.6). Thus, 



(37) ^+=0+^/ dr, 



, no . nu> a f°° , 1 1 a 2 f°° , 1 1 1 

(38) Q = 2tt ]_„ ^W^WT* + ^WT^D^TV^WTVr, 

where the first summand of the right hand side is denoted by aQ\ and the second 
summand is treated in ljl4|l through (|16|> . 
Furthermore, since 

d P ^mi <00 
i+ipi 

the potential ip is regular in the sense of Klaus and Scharf |15| , namely the operator 
Q v £ & z (Sj), i.e., Q v is an Hilbert-Schmidt operator. (See also [T2*]. Theorem 
4.) This allows to show the finiteness of the exchange energy between the one- 
particle density matrix of the electron-positron field and the difference of the density 
matrices of the polarized Dirac sea and the free Dirac sea. 

To formulate the next lemma we fix the following notation: let C v%q be the 
optimal constant in the generalized Young inequality, i.e., ||/*<?||r < Cp.gll/llpllsllgj 
1 < p,q,r < 00, r _1 + 1 = p^ 1 + q^ 1 . 
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Lemma 2. Let tp e L 3 (T) n L 2 (T). Then 
(39) 



dx / dy 



rl>(x)Q*(x,y)1>(v) 



< A/ C3/2,3/2ll 1/ 



and /or every e > i/iere exists a constant C e > smc/i £/iai 



13/2,- 



(40) 



/ dy 



il>{x)Q*{x,y)^{y) 



<eU\\l + C e U\£ 



l*-y| 

Proof. Since Q v is a Hilbert-Schmidt operator we get using the Schwarz inequality 



(41) L:= 



dx / dy 



|x-y 



< ( dx dy 



l x -y| : 



2\ 1/2 



dx y dy|Q ¥> (x,;. 



1/2 



The second factor of the right hand side is the Hilbert-Schmidt norm || || 2 of Q* '■ 
To estimate the first factor we decompose the kernel into two functions /(x) := 

XBr(o)( x )/I x | 2 and the rest g, i.e., l/|x| 2 = /(x) + .g(x). 
Thus, using inequality (|41|l we get 



(42) 



L < 



/,M 2 ) 1/2 + (H 2 *s,H 2 ) 1/2 



We estimate the first and second summand of the first factor on the right hand side 
separately from above. 

The first summand yields using the Holder inequality followed by the generalized 
Young inequality (see, e.g., Reed and Simon [22], P- 32) 



(43) 



/,|^| 2 )<C 3A3/2 ||^ 2 ||2 /2 ||/||3 /2 ^, 



where w indicates the weak-norm. 

Picking the radius R — 00, i.e., 5 = 0, yields immediately (|39|l . 

To prove (|40|l we also use (|43|l but pick the radius R > sufficiently small: in 
this case we need to bound also the second summand containing <?; we use again 
Holder's inequality now followed by using Young's inequality 

(44) (IVf ^H 2 )^^ 2 !! 2 !!^- 

Thus, the first factor on the right hand side of l|42|l is bounded by 



3/2,3/21 



ll^ll3/2ll/ll3$L + IIV- 2 ||l||3llL /2 - 



Since II/H3/2,™ tends to zero as R tends to zero, the claimed inequality follows. □ 



3.2. Electrostatic Vacuum Polarization Energy (Non-Exchange Energy). 

In the expression for the electrostatic vacuum polarization energy we replace the 
density of the polarized sea by the difference of this density and the free Dirac sea 
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using Principle W2: 

J J |x-y| 
(45) =Jdxjdy: **(a?)*(x) : T^jr 

(Here and in the following we will denote by Q v (x, y) the 4x4 matrix with entries 
(Q v (x, er; y, r)) r _ r ) However, the integral kernel of Q v is always singular on the 
diagonal except for vanishing potential as can be seen from i|51[l implying that P 
is not well defined; one more renormalization is necessary. The question how to 
extract the physical relevant information from ti^Q^y, y) was already asked by 
Dirac [2] and partially answered by Dirac [3], Heisenberg [TT], Serber [ZB|. Uehling 
|2"§| . Weisskopf [3T], Schwinger |2S], Dyson 0], Klaus and Scharf an d others. 
The proposed solution amounted to a perturbative renormalization according to 
Principle W3. — One of our main results is that this renormalization can be done 
non-perturbatively: subtracting the zeroth order expansion Q\ of the difference Q v 
of and P® will turn treaty (y, y) into well defined quantities given in (|f 9JI and 

ipnji. 

Recall that 

1 f°° 1 1 



2tt7_ 00 'D° + irj D° +ir)' 

Thus, in momentum space Q\ is given by 

a ■ p + j3 — irj ^ _ N a-q + /3 — i?7 
q 2 + 1 + ?/ 2 



(47) QrC^q) = (2^)" s / 2 f - P + % (p - q) 



which leads to 

(48) trc&Cp, q) = 2-V^-»/V(p - q) P; q +;-f ( Pf^\, 



by a straightforward calculation with E(p) = -\/p 2 + 1. In configuration space we 
obtain 

(49) 

tr C 4Q 1 (x,y) = (2 7 r)- 3 / dr / dqe <Mt tr C 4&(r ) q) e - i «« r 

= (2tt)- 3 /" dpdktr C 4Qi(p-k/2,p + k/2)e <p -( x - y )e- ik '( x+y) / 2 
=:Q(x- y ,^) 

after introducing new variables of integration r = p— k/2 and q = p+k/2. Defining 
(:=x-y we remark that the "limits" lim y ^ x trc4<2i(x, y) and lim^o Q(£, x) are 
formally the same. The corresponding expression P in the electrostatic energy Q45JI 
becomes formally 

(50) J ix J Jy - «• W W r) _ J w g ftp, k) , 
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where Q(£, •) is the Fourier transform of Q with respect to the second variable for 
fixed £ ^ 0, i.e., formally 

(51) Q($, k) = (2tt)- 3 / 2 j ^ dp tr C 4Q 1 (p - k/2, p + k/2)e lp « 

= J-^fk) / rfp P 2 - k2 / 4 + 1 - g(E - k / 2 )^P + k / 2 ) c ** 
4tH» W V R 3 P £(p-k/2)£(p + k/2)(£(p-k/2) + i;(p + k/2)) 

We note that the integral l|51l) is logarithmically divergent at £ = independently 
of the form of the external potential tp. This shows - as already remarked above - 
that the limit lim y _> x trc<iQ v (x, y) only exists, if <p vanishes. 

The expression Q(£,k) has been intensively studied in the literature (see, e.g., 
Heisenberg Serber (2E|, Pauli and Rose |2()j . Weisskopf (HJi an d Klaus and 
Scharf 16 ). We will follow mainly the calculations of Pauli and Rose. (Since their 
treatment is time-dependent one has to set ko = to translate to our situation.) 

According to [301, Equations (5) - (9), we can separate Q into two terms 

(52) ^(^k)=^ 1 (^k) + ^(k)k 2 J Fb(0. 
with 

1 f , ( p 2 cos 2 6>\ e ip < 

(53) F (£)=_ / dp 



16tt 3 J m z\ 1 + p 2 J (1 + p 2 ) 3 / 2 

9 being the angle between £ and p. With this definition of Fq the function F% is 
finite for £ = and has there the value 



(54) p vac (k) :=fi(0 t k) 

_y(k) f p 2 -^ + l-g(p-§)£(p + |) k2 p 2 sin 2 g + l 

4^ 3 7 K 3 E(p - |)£(p + f ) (E(p -|) + E{p + I)) 4£(p)5 



rfp 



where C is the function defined in 1)21(1. While each of the summands in the latter 
formula decreases like | p> | 3 for large values of |p| and therefore the corresponding- 
parts of the integral are logarithmically divergent, the difference in the integrand 
decreases like |p| -5 and is therefore convergent. 

Pauli and Rose [20] obtain the following asymptotic behavior for C: 



(55) C(k)/k 2 = 



2 J 15 



i,k 2 +o(k 2 ) | k | 



2 log(|k|)-|+ (l) |k| 



We note that the second summand on the right hand side of (|52|l can be written 
as 47rn(k)i<o(£), i.e., this depends only on the density of the nucleus. This implies 
that it can be dropped according to W3. This means that p vac as defined in lf54T) 
can be considered as the physically relevant density of the polarized vacuum. 

3.3. The Fully Renormalized Hamiltonian. The potential of the vacuum po- 
larization density is 



(56) U=\-\- L *p- 



vac i 
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the corresponding term of the electrostatic interaction of the electron-positron field 
with the vacuum reads 



(57) -cr/ dx : **(x)*(a:) : U(x) = -a 1 J dx : : | • \~ L *p vac (x) 

= -a 2 Jdk:^* : (fc)^ Pvac (k) 

where we used .F(| • | _1 ) = y/2/n\ ■ |~ 2 . 

According to (|45() , (|38|l , l|17|l , and 1|19|) the second quantized renormalized polar- 
ization energy becomes 

(58) P = J dx: : (f7(x) + a 2 P 3 (x) + a 3 F 4 (x)) . 
Consequently our fully renormalized Hamiltonian is 



(59) I = F + aW + a 2 (-P + X) 



3.4. Physical Interpretation of the Renormalization Procedure of the 
Vacuum Polarization. In physics literature the subtraction of the singular part 
of the diagonal term of the one-particle density matrix of the Dirac sea, i.e., the 
dropping of the second summand ]b'2\i ) is called charge renormalization for the 
following reason: the term subtracted from P to obtain P is 

(60) a 2 47rF (£) J dx : V*(x)V(x) : <p(x). 

Formally this result can also be obtained by replacing the square of the (bare) 
charge e 2 = a in D v 

(61) e 2 i-> e 2 (l -ATre 2 F {£)). 
Note, that 

F (£) = -log(|£|) + O (l) 

for small |£| which leads to a well known formula in the literature (see, e.g., Milonni 
|19j . p. 417). It is interesting that a change in the effective charge due to the 
polarization of the vacuum was already suggested by Furry and Oppenheimer jS] . 

The word polarization is due to the following picture: according to Dirac the 
electrostatic field causes a redistribution of charge in the Dirac sea, i.e., it polarizes 
the vacuum. In particular the nucleus polarizes the vacuum in its vicinity causing 
a screening lowering the effective charge for an observer at a distance. 

Another reason for the fact that the infinity of the diagonal part of the density 
matrix of the Dirac sea in (|52(l is invisible in experiments is the following: in first 
order the factor in front of Fq(0) changes a in D v by an (infinite) constant only, 
which does not effect the degeneracy of the eigenvalues of D v , i.e., it does not cause 
a splitting of degenerated eigenvalues, a fact that is confirmed by experiment. 

3.5. The Vacuum Polarization Potential of the Coulomb Potential. Recall 
that the nuclear potential is ip — \ ■ | _1 * n; consequently ip(k) = 47rn(k)/k 2 . Thus 
the Fourier transform of the vacuum polarization potential using (|54|) gives 

ttm ^W c ( k ) X k )g( k ) 
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which is spherically symmetric and compactly supported under our assumptions on 
the charge distributions of the nucleus. 

If we assume that the nuclear density is a spherically symmetric Schwartz func- 
tion - as we do in the mathematical part of this paper - this implies that the vac- 
uum polarization potential U is bounded continuous and decreases exponentially 
at infinity. 

However, if we assume that we have a point nucleus as assumed by Uehling |29| . 
this is no longer the case. To relate to Uehling's work we will discuss this case as 
well although the corresponding potential is no longer form bounded with respect 
to the kinetic energy. 

From l)62|) we have 

(63) U = Z^T-^ 

According to Uehling (2H| and Schwinger Equation (2.53), this is 

which means asymptotically 

to* r7 ^_/-^W- 1 0°8W + i + 7) + o(i) M-o 

( j { ) ~ l^e-'Mlxl^Cl + 0(l/|x|)) |x| - oo 

where 7 w 0.5772 is Euler's constant. Consequently one obtains an effective poten- 
tial 

Z 2 Z f 5 \ 

(66) ^eff(x) = -cap - a 2 U = -a— + a 2 — — log |x| + - + 7 + 0(1) 

|x| 3tt |x| V 6 7 

close to the nucleus. 

Obviously, due to vacuum polarization, the effective potential becomes more 
singular than the Coulomb potential. This implies that the energy 

II 2 2Z 1 
3tt|x| |x| 

is unbounded from below for all positive values of a and Z and l|66|) is no longer 
relatively form bounded with respect to the relativistic kinetic energy operator. 
This suggests to avoid the mathematical idealization of a point nucleus and to take 
the experimental fact that the nuclei are extended into account. 

3.6. Splitting of the Bound State Energies. The effect of the vacuum polar- 
ization potential to lowest order in a is given by the effective one-particle operator 

P|(^-a 2 f/(x))^. 

Therefore, the energy eigenvalues are shifted in lowest order of a by 

(67) 8E = ~a 2 I dx[/(x)|^(x)| 2 

where ip(x) denotes an eigenstate of D v . 

To get a rough heuristic estimate on the numerical effect Uehling assumes 
the nucleus to be a point particle, i.e., its density is n(x) = Z5(x), and takes 
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the corresponding Schrodinger eigenstates ip n ,i where n is the principal quantum 
number and / the orbital- angular- momentum quantum number: 



(68) 5E n j = -a 1 / dx[/(x)|^,z(x)| 
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l^n,;(0)| 2 = - 



4^,5 



4Z 4 a 
157rn ; 



Concerning the first excited eigenvalue n = 2 this indicates an energy level 
splitting (of the 2s and 2p state) of 



(69) 



53 



2,0 



Z^oc'ra 
30 ' 



the Uehling effect [211, Equation (29). 

The vacuum polarization (|69J) accounts for only one percent of the 2si/2 — 2pi/2 
Lamb shift of hydrogen, since the Bohr radius is much bigger than the range of 
the vacuum polarization potential. However, the Bohr radii of muonic atoms are 
much smaller because of the large effective mass which means that the vacuum 
polarization of muonic helium accounts for 90 percent of the Lamb shift (Peterman 
and Yamaguchi [2*T]. Glauber et al. [5], see also Greiner et al. JU|j P- 413). 



4. Proof of the Self-Adjointness of the Renormalized Hamiltonian 

In the following we are going to show that the quadratic form associated to the 
operator H of the electron-positron field restricted to the one-electron sector $j+ 
defines a self-adjoint operator which is bounded from below. This means - among 
other things - that higher order renormalizations, as introduced in perturbation 
theory by Dyson 0j, are unnecessary. 

We will show infinitesimal form boundedness of all perturbations, namely U, P3, 
and P4 relative to P^D^P^. For X this has been already shown in Lemma|21 

U: According to H62|) and the remark thereafter U is bounded and therefore 
trivially infinitcsimally relatively form bounded. 

P4: We will show that for any positive e there exists a C e such that for all 

ip e v(d^) c L 3 (r) nL 2 (r) 



(70) 



{^P^)<eU\\l + C t U\\ 



which implies the infinitesimal form boundedness by the Sobolev's inequality for 



Lemma 3 

(71) 
with 
(72) 

Proof. We have 



Assume aZ 6 [0, 1), x € ^ 5 (r)- Then 



tr(xQ 4 )| < \\xQ4 



< 



3tt 2 



<P\\5\\X\\5 



a 



(73) 



tr( X Q 4 )| <5- 



<- 



'2vr 



dr/ 
dr/ 



1 



H-«ll^l 



1 



1 



-f- 



D° + ir) r D 
1 



D° +IT] 



— w 

■ir) Dt> 



D° +irj 



-V 



IT] 



D° 
1 



IT) 



D° 



IT] 
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(We use the standard notation ||A|| P = ^/tr \ A\p.) We will estimate the right hand 
side of the above inequality which will also show that xQ^ £ 61(f)). To this end 
we estimate the factor containing the perturbed resolvent: 

(74) y-J—\\ 5 = y-^—(D° + z V )- ' 



D'P + ir]" " D° + ir) K u Df + i-q 



< ||(J> +^) n[O , J Uly- 



The first factor on the right side is finite: 

UD ° + tv) D^Tu^ ~ l + a ^D^ ] U < l + «ll^lloo <oo 

where we use that D v is invertible because of Lemmas and that ip is bounded (see 
Q). (Note the boundedness would also hold if <p = Za/\ ■ |, since l/ \ ■ | is relatively 
bounded with respect to >/—A.) 
Since 

(75) \\f(x)g(- l \7)\\ 5 <(2nr^\\f\\ 5 \\g\\ 5 , 

the norm being the one of the trace ideal 6s(L 2 (M 3 )) (Simon [57], Theorem 4.1), 
we can estimate the other factors occurring in J73J): 



2 



5 



(76) "X5o+^ll« ^ W 5^llxll 5 ||l/V|-| 2 + 1 + ^ 

(norm in ©5 (55) on the left hand side and in S5(L 2 (R 3 )) on the right hand side) 
and a similar expression of the term containing (p. 
Using (|74*jl and (f7H|l allows us to continue lf?3"|) as 

n r°° r 1 

11 11 11 ||4 / j I j 



f t />oo p 



(p 2 + l+T 1 2 f/ 2 

CipA f j 1 f 1 , I! I! I! 11 4 ^fA II 11 11 ||4 



□ 



This has the 



Corollary 2. TTie perturbation P4 is relatively form bounded with respect to \D°\ 
with form bound zero. 

Proof. We pick x = l^l 2 * I ' I 1 m Lemma |3] with i/> £ i? 1 / 2 (r). Using Young's 
inequality followed by Sobolev's inequality yields the desired result. □ 

P3: Unfortunately, Simon's elegant trace inequality used in (|73|l does not suffice 
to handle the Q3 containing only four resolvents. In that case we estimate directly: 

Lemma 4. Denote the electric potential of a state ip S H 1 / 2 ^) of finite kinetic 
energy by x( x ) := f dy\ip(y)\ 2 /\*- — y| and assume p > 3. Then there exists a 
constant C' v-P such that 

(77) \(M^)\<o v Ax\\p- 
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Proof. We have 
(78) ty,P#l>)= / dx X (x) P3 (x) 



dPX(p)P3(p) 



(2tt) 3/2 / dpi / dp 2 Vx(Pi ~ P2)Q3(Pi,er;P2,cr) 

JR3 J R 3 ct=1 



where we use the Definition (|17|l of p 3 . The "eigenfunctions" of the free Dirac 
operator in momentum space are 



(79) 



u T (p) := < 



Tf+W 



a ■ pe T 

-(l-£(p))e T/ 
cr • pe T 



r= 1,2, 



r = 3,4 



\-(l + E(p))e T/ 
with e T := (1,0)* for r = 1, 3 and e T := (0, 1)* for r = 2, 4 and 
(80) JV+(p) = V2^(p)(£;(p) - 1), N-(p) = y/2E(p)(E(p) + l). 

The indices 1 and 2 refer to positive "eigenvalue" B(p) and the indices 3 and 4 to 
negative — E(p), (See, e.g., Evans et al. UJ.) Using Plancherel's theorem we get 



(81) (V> 



P 3^) = TTpry / dpi / dp 2 [ dp 3 [ dp 4 V 

v / ^R 3 JR 3 JR 3 JR 3 ri!T3 ^ T4=1 



X(Pl - P2)<^(P2 - P3)^(P3 - P4)<^(P4 - Pl)x 
X K (Pl)|«T 2 (P2)) ( u r 2 (P2)|«T 3 (P3)) («t 3 (P3)|"r 4 (p<t)) («r 4 (P4)|«n (Pi)) 

1 



di] 



(ia Tl E{pi) - r))(ia T2 E(p 2 ) - ri)(ia T3 E{p 3 ) - rj){ia Ti E(p 4 ) - i]) ' 



with a T — 1 for r = 1, 2 and a T = — 1 for r = 3, 4. The integral over r/ is seen to 
vanish by Cauchy's theorem, if all four a Tj have the same sign. In fact we have to 
distinguish only two cases, namely three of the a T . are equal and two of the a T - are 
equal. 

Therefore in the following we will only treat two different cases. The others then 
work analogously. 
We begin with 



(82) a Tl 
In that case the first factor in (|81|l reads 

(83) K 1 (Pl)|Wr 2 (p 2 )) K 2 (P 2 )K 3 (P3)>X 

7-1=3,4 7-2 = 1,2 

X Y ( U 7- 3 (P3)K 4 (P 4 )) Y K 4 (P4)Ki(Pl)) = 

7-3 = 1,2 7-4=1,2 

cr • pier • p 2 + (1 + £(pi))(l - E(p 2 )) 



tr c 



■ [a • p 2 G • P3 + (1 - E{p 2 ))(\ - E(p 3 ))] 



iV_(p 1 )2A/ + (p 2 )2A r+ (p 3 )2 Ar+( p 4) 2 

[a • p 3 c7 ■ p 4 + (1 - E(p 3 ))(l - E( Pi ))] [a ■ p 4 cr • pi + (1 - E( Pi ))(l + E( Pl ))] 
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We estimate the modulus of (|83J) and obtain 

tr C 2 \a ■ p 4 a • pi + (1 - E(p 4 ))(l + £(pi))| 



IdHHH < c- 

< c- 



iV_( Pl )iV + (p 4 ) 
p 4 Pi - (E(p 4 ) - 1)(1 + £( Pl ))| + |p 4 A Pl | 



7V_( Pl )7V + (p 4 ) 

(Here and in the following c is a generic positive constant.) Since 

j_ r d i 

1 



(E(p 2 ) + E{p 3 ))(E(p 2 ) + E(p 3 ))(E(p 3 ) + E(p 4 )) 
our term of interest l|81|) is bounded by constant times 

(84) 

/ dpi / dp 2 dp 3 dp 4 |x(Pi ~ P2)^(P2 - P 3 )<£(p 3 ~ P4)^(P4 - Pi)| 

Jr 3 Jr 3 Jr 3 Jr 3 

|p 4 • pi - (E(p 4 ) - l)(£(pi) + 1)1 + |p 4 A pil 



N^(p 1 )N+(p 4 )(E(p 2 ) + E(p 3 ))(E(p 2 ) + E(p 3 ))(E(p 3 ) + E(p 4 ) 
Substituting p 2 — * Pi + p 2 turns into 

(85) / dpi / dp 2 / dp 3 / dp 4 

JR 3 JR 3 JR 3 Jr 3 

\if{p 2 + Pi P 3 )<^(P 3 - P4,)<f(Pi ~ Pl)| 



)■ 



Ix(-p 2 )l 



N-( Pl )N+(p 4 ) 
|p 4 ■ Pi - (E(p 4 ) - l)(g(pi) + 1)| + |p 4 A pi| 
{E(p 2 + pi) + E(p 3 ))(E(p 2 + Pl ) + E(p 3 ))(E(p 3 ) + E( Pi ) 

^P2|X(-P2)|/(P2) 



where we introduce / to be the remaining integrand. We will now estimate /. 
Substituting pi — > pi + p 4 , p 3 — > p 3 + p 4 we get 

/(P2) = / dpi / dp 3 dp 4 <p(p 2 + pi - p 3 )^(p 3 )^(pi)| 

JR 3 JR 3 JR 3 

x |P4 • (Pi + p 4 ) - (E(p 4 ) - 1)(1 + Ejpi + p 4 ))| + |p 4 A pil 
iV_(pi + p 4 )iV + (p 4 )(£:(p 2 + pi + p 4 ) + E(p 3 + p 4 )) 

1 



(E(p 2 + pi + p 4 ) + E(p 3 + p 4 ))(E(p 3 + Pi ) + E(p 4 )) ' 

Since 

S(Pl+P4)=^(P4) + ^7-^yM 

for some /j, £ [0, 1], we see that 

|P4 • (Pi + p 4 ) - (E(p 4 ) - l)(E(p 1 + p 4 ) + 1)| + |p 4 A pi| < 4|pi||p 4 |. 
Notice, we can bound 



/ 

Jr 3 



, |P4| 

ap4 E(p 3 + p 4 )E(p 4 yN + ( Pl + p 4 ) - c 
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independent of pi and P3. Therefore, 

/(P2) < c / dpi / dp 3 |0(pa +Pl -P3)0(P3)<£(Pl)||Pl|- 

Since </?(k) = 47m(k)/k 2 , we have that /(0) is finite; since n is compactly supported, 
Pi and p 2 are bounded. We conclude that / has also compact support. 
Consequently, since 

(86) I dp 2 |x(-p 2 )|/(p 2 ) < ||xll,ll/llp, 

we see using the Hausdorff- Young inequality that for all p > 2 



(87) 



/ dp 2 |x(-p 2 )|/(p 2 ) < c PtV \\x\ 
Jr 3 



for constant c p . v depending on p and (p. 

Next, we take a peek at the case a Tl = a T2 = 1 and a T3 = a T4 = — 1. The 
corresponding integral over n gives 



2nJ_ 00 (iE( Pl ) - vKiEfa) - ri)(-iE(ps) - fj)(-i£7(p 4 ) - v) 

1 

(E(p 2 ) + E(p 3 ))(E(p 2 ) + E(p 4 ))(E( Pl ) + E(p A j) 

1 



(E(p 2 ) + E(p 3 ))(E( Pl ) + E(p 3 ))(E( Pl ) + E(p 4 )) 

Observe now that the corresponding first factor in l)81|l can be bounded by c ■ 
4 AT J P2 il^ 3 l — 7 • Now, we do similar variable transforms as above and arrive at an 

JV + (p2)A f -(p 3) ' 

analogue of (|57Jl . □ 
Again this has the 

Corollary 3. The perturbation P 3 is relatively form bounded with respect to \D°\ 
with form bound zero. 

The above result was the final step in showing that all four perturbations, U, 
P3, P4, and X, are form bounded with respect to PfD v Pf. 



References 

[1] H.A. Bethe. The electromagnetic shift of energy levels. Phys. Rev., Minneapolis, II. Ser., 
72:339-341, 1947. 

[2] P.-A.-M. Dirac. Theorie du positron. In Cockcroft, J. Chadwick, F. Joliot, J. Joliot, N. Bohr, 
G. Gamov, P.A.M. Dirac, and W. Heisenberg, editors, Structure et proprietes des noyaux 
atomiques. Rapports et discussions du septieme conseil de physique tenu a Bruxelles du 22 
au 29 octobre 1933 sous les auspices de I'institut international de physique Solvay. Publies 
par la commission administrative de I'institut., pages 203-212. Paris: Gauthier-Villars. XXV, 
353 S., 1934. 

[3] P.A.M. Dirac. Discussion of the infinite distribution of electrons in the theory of the positron. 

Proc. Camb. Philos. Soc, 30:150-163, 1934. 
[4] F. J. Dyson. The radiation theories of Tomonaga, Schwinger, and Feynman. Physical Rev. 

(2), 75:486-502, 1949. 

[5] F. J. Dyson. Divergence of perturbation theory in quantum electrodynamics. Physical Rev. 
(2), 85:631-632, 1952. 



VACUUM POLARIZATION 



19 



[6] William Desmond Evans, Peter Perry, and Heinz Siedentop. The spectrum of rclativistic one- 
electron atoms according to Bcthc and Salpctcr. Comm. Math. Phys., 178(3):733-746, July 
1996. 

[7] J.D. French and V.F. Weisskopf. The Electromagnetic Shift of Energy Levels. Phys. Rev., II. 

Ser., 75:1240-1248, 1949. 
[8] W.H. Furry and J.R. Oppcnhcimcr. On the Theory of the Electron and Positive. Phys. Rev., 

II. Ser., 45:245-262, 1934. 
[9] Roy Glauber, William Rarita, and Philip Schwcd. Vacuum polarization effects on energy 

levels in /i-mesonic atoms. Physical Review, 120(2):609-613, October 1960. 
[10] W. Greiner, B. Miiller, and J. Rafelski. Quantum Electrodynamics of Strong Fields. Texts 

and Mongraphs in Physics. Springer- Verlag, Berin and Heidelberg and New York and Tokyo, 

1 edition, 1985. 

[11] W. Heisenberg. Bemerkungen zur Diracschen Theorie des Positrons. Z. Phys., 90:209-231, 
1934. 

[12] Bernard Hclffcr and Heinz Siedentop. Form perturbations of the second quantized Dirac field. 

Math. Phys. Electron. J., 4:Paper 4, 16 pp. (electronic), 1998. 
[13] J. M. Jauch and F. Rohrlich. The theory of photons and electrons. The relativistic quantum 

field theory of charged particles with spin one-half. Addison- Wesley Publishing Company, 

Inc., Cambridge, Massachusetts, 1955. 
[14] Tosio Kato. Perturbation Theory for Linear Operators, volume 132 of Grundlehren der math- 

ematischen Wissenschaften. Springer- Verlag, Berlin, 1 edition, 1966. 
[15] M. Klaus and G. Scharf. The regular external field problem in quantum electrodynamics. 

Helv. Phys. Acta, 50(6):779-802, 1977. 
[16] M. Klaus and G. Scharf. Vacuum polarization in Fock space. Helv. Phys. Acta, 50(6):803-814, 

1977. 

[17] Norman M. Kroll and Willis E. Lamb jun. On the Self-Energy of a Bound Electron. Phys. 
Rev., II. Ser., 75:388-398, 1949. 

[18] Willis E. Lamb and Robert C. Retherford. Fine structure of the hydrogen atom by a mi- 
crowave method. Physical Review, 72(3):241— 243, August 1947. 

[19] Peter W. Milonni. The Quantum Vacuum: An Introduction to Quantum Electrodynamics. 
Academic Press, Inc., Boston, 1 edition, 1994. 

[20] W. Pauli and M.E. Rose. Remarks on the Polarization Effects in the Positron Theory. Phys. 
Rev., II. Ser., 49:462-465, 1936. 

[21] A. Peterman and Y. Yamaguchi. Corrections to the 3d-2p transitions in /i-mesonic phosphorus 
and the mass of the muon. Physical Review Letters, 2(8):359-361, April 1959. 

[22] Michael Reed and Barry Simon. Methods of modern mathematical physics. II. Fourier anal- 
ysis, self-adjointness. Academic Press [Harcourt Brace Jovanovich Publishers], New York, 
1975. 

[23] Michael Reed and Barry Simon. Methods of Modern Mathematical Physics, volume 4: Anal- 
ysis of Operators. Academic Press, New York, 1 edition, 1978. 

[24] Silvan S. Schwcbcr. QED and the men who made it: Dyson, Feynman, Schwinger, and 
Tomonaga. Princeton Series in Physics. Princeton University Press, Princeton, NJ, 1994. 

[25] Julian Schwinger. Quantum Electrodynamics II. Vacuum Polarization and Self-Energy. Phys. 
Rev., II. Ser., 75:651-679, 1949. 

[26] Robert Serber. Linear modifications in the Maxwell field equations. Phys. Rev., II. Ser., 
48:49-54, 1935. 

[27] Barry Simon. Trace Ideals and their Applications, volume 35 of London Mathematical Society 

Lecture Note Series. Cambridge University Press, Cambridge, 1979. 
[28] Bcrnd Thaller. The Dirac Equation. Texts and Monographs in Physics. Springer- Verlag, 

Berlin, 1 edition, 1992. 

[29] E.A. Uehling. Polarization effects in the positron theory. Phys. Rev., II. Ser., 48:55—63, 1935. 

[30] Steven Weinberg. The quantum theory of fields. Vol. I. Cambridge University Press, Cam- 
bridge, 1996. Foundations, Corrected reprint of the 1995 original. 

[31] V. Weisskopf. Uber die Elektrodynamik des Vakuums auf Grund der Quantcntheorie des 
Elektrons. Math.-Fys. Medd., Danske Vid. Selsk., 16(6):l-39, 1936. 



20 C. HAINZL AND H. SIEDENTOP 

Mathematik, Theresienstrasse 39, 80333 Munchen, Germany. Address as of March 
1, 2003: CEREMADE, Universite Paris- Dauphine, Place du Marechal de Lattre de Tas- 
signy, F-75775 Paris Cedex 16, France. 

E-mail address: hainzlSceremade.dauphine.fr 



Mathematik, Theresienstrasse 39, 80333 Munchen, Germany. 
E-mail address: h. s31mu.de 



